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WITHOUT A RADIAL ASSUMPTION 
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Abstract. We extend our result [29] to the non-radial case, giving a complete 
classification of global dynamics of all solutions with energy at most slightly above 
that of the ground state for the nonlinear Klein-Gordon equation with the focusing 
cubic nonlinearity in three space dimensions. 



1. Introduction 

The nonlinear Klein-Gordon equation 

Uu + u = ±u 3 , (t, x) e M.]+ 3 (1.1) 

has been considered by many authors. The signs ±u 3 are known as focusing (+) 
and defocusing (— ), respectively. This reflects itself in the conserved energy 



E(u, u) 



^(|Vn| 2 + \u\ 2 + |n| 2 )^u 4 dx. (1.2) 



In either of these cases, short-time existence of smooth solutions is known, as well 
as global existence and scattering to zero for small data. For the defocusing case, 
one has global existence for all data as well as scattering in the energy class 

u(t) := (u(t),u(t)) G H := H 1 x L 2 , (1.3) 

see Brenner 0HD], Ginibre, Velo [161 [T7], Morawetz, Strauss [28], and Pecher [32] . 
As usual, scattering of a solution u to a static state cp refers to the following asymp- 
totic behavior: there exists a solution v of the free Klein-Gordon equation in the 
energy class such that 

\\u(t) - <p - v(t)\\ n ^ 0, (1.4) 

as either t — > oo or t — > — oo (depending on the context). When ip = 0, we sometimes 
say "scattering" instead of "scattering to 0". See also Strauss [36] and [19] for a 
review of Strichartz estimates and wellposedness, as well as scattering in this setting. 
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dx 



In this paper, we are only concerned with the focusing equation, for which smooth 
solutions can break down in finite time (which means that the L°°-norm becomes 
infinite in finite time). More generally, the blow-up is defined for weak solutions in 
the energy class by divergence of the energy norm % (which is more natural in view 
of the wellposedness theory, but these two definitions are equivalent for classical 
bounded^ solutions). In fact, this happens for any data of negative energy, as was 
shown by Levine [25] . A natural question is then to decide which data lead to finite 
time blow-up versus scattering to zero. A first step in this direction was undertaken 
by Payne and Sattinger [31]. Their result is formulated in terms of the so-called 
ground state, which is the unique positive radial decaying solution of 

-AQ + Q = Q 3 , 

see [351 [?l [11]. Amongst all stationary solutions of (II. ip in H 1 , ±Q are unique (up to 
translations) with the property that they minimize the action (or stationary energy) 

A scaling functional associated with J is defined by 

k(<p) = [ V ^I 2 + M 2 - Ml dx 

The significance of K lies with the fact that 

inf{ J(<p) | K(<p) = 0, y? e H 1 \ {0}} = J(Q) 

As observed by Payne and Sattinger, this easily implies that the regions 
VS + = {{u , Ul ) E H | E(u , Ul ) < J(Q), K(uo) > 0} 
VS^ = {(m , mi) G H | E{u , Ul ) < J(Q), K(u ) < 0} 

are invariant under the nonlinear flow in the phase space %. Moreover, they showed 
that solutions in VS + are global (in both time directions), whereas those in VS_ 
blow up in finite time (in both time directions). In particular, the stationary solution 
Q is unstable, see also Shatah [31] and Berestycki, Cazenave [6]. Scattering in VS + 
was only recently shown by Ibrahim, Masmoudi, and Nakanishi [19], using a Kenig- 
Merle type argument originating from [221 123] . 

This paper addresses the behavior of solutions whose energy exceeds that of J(Q) 
by a small amount. More precisely, let u := (u, u) and set 

U e := {u e U | E m (u) < E(Q,0) + e 2 } (1.6) 

with the minimal energy 

E m (u) := \E(u) 2 - P{uf\^ sign(£(n) 2 - P(n) 2 ), (1.7) 



1 One can easily construct smooth global solutions with finite energy norm but with infinite L° 
norm from the beginning, by superposing very sparse sequence of concentrated waves. 
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where P(u) = (ii|Vu) is the conserved momentum. In contrast to E, the minimal 
energy E m is Lorentz invariant, which makes it the right quantity to use in the 
nonradial case. We call a solution with zero momentum normalized. To every 
solution with |P(w)| > |-P(w)| there exists exactly one Lorentz transform^ which 
reduces it to a normalized one, see (I3.2p . If l-E^u)! < |P(w)|, then there exists a 
Lorentz transform L such that E(uoL) < 0. But such solutions are known to blowup 
in finite positive and negative times, see Payne-Sattinger [31]. If | -£?(«)[ = |P(u)|, 
then along some sequence of Lorentz transforms Lj we have E(uoLj) — > as j — > oo. 
But then either K(u o Lj) < for some j, which means that u blows up in both 
time directions by (21] , or u o Lj — > strongly in H 1 so that u exists globally. Thus, 
we can always talk about normalized solutions in those cases where [3T] does not 
apply, which are the ones relevant to this paper. 

Theorem 1.1. Consider all solutions of NLKG fll.ip with initial data u(0) e 1-L e for 
some small e > 0. We prove that the solution set is decomposed into nine non-empty 
sets characterized as 

(1) Scattering to for both t — > ±oo, 

(2) Finite time blow-up on both sides ±t > 0, 

(3) Scattering to as t — >■ oo and finite time blow-up in t < 0, 

(4) Finite time blow-up in t > and scattering to as t — > — oo ; 

(5) Trapped by ±Q for t — > oo and scattering to as t — > —oo, 

(6) Scattering to as t — >■ oo and trapped by ±Q as t — > — oo, 

(7) Trapped by ±Q for t — > oo and finite time blow-up in t < 0, 

(8) Finite time blow-up in t > and trapped by ±Q as t — >■ — oo, 

(9) Trapped by ±Q as t — > ±oo, 

where "trapped by ±Q" means that the normalized solution stays in a 0(e) neigh- 
borhood of 

{±Q(- + y)\yeR 3 } 

forever after some time (or before some time). The initial data sets for (l)-(4), 
respectively, are open. 

The radial case of this theorem was proved in [29] . We follow the same approach 
as in [29], but with the added feature of having to control a small translation vec- 
tor, after the Lorentz transform to the normalized solution. As in the radial case, 
Theorem 11.11 extends to all dimensions x G M. d and equations of the form 

4 4 
Du + u = |w| p "V l + ^<p<l + ^ — -, 

a a — 2 

namely those powers of L 2 -supercritical and iP-subcritical. The main difference 
between [29] and this one lies with the further analysis of the trapped solutions 
in Theorem 11.11 Recall that we proved in the radial case (where all solutions are 



For subtle issues such as the meaning of the Cauchy problem for Lorentz transformed solutions, 
and general wellposcdncss questions of transformed solutions, see Section [S] below. 
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normalized) that any solution in T-L e which is trapped by Q as t — > oo scatters to Q 
as t — > oo, see (11.41) . 

This is the "conditional asymptotic stability behavior" observed in [551 I2"4"l [5] for 
the nonlinear Schrodinger equation (NLS), and near the equilibrium Q these states 
form a smooth codimension one manifold which is of the center-stable type, see [3J . 
Note that there are no modulation parameters in the radial case, such as translation 
or scaling, which simplifies the analysis considerably. In contrast, in the nonradial 
case, Q may be translated as well as Lorentz transformed, leading to a 6-dimensional 
manifold of solitons (which consists of two disconnected components corresponding 
to ±Q), parametrized by the relativistic momentum peR 3 and the center of mass 
q G M 3 . The traveling waves also undergo the Lorentz contraction (i.e., they are 
flattened in one direction), and are explicitly given at each (p, q) G M 6 in the form 

Q(p,q) = Q(x -q + p((p) -i)bl"V p>^i 3 , (1.8) 

where (p) : = \/\ + \p\ 2 . The ground-state traveling waves are those solutions 

u{t) = ±Q(p, q(t)), p G R 3 , q{t) = ^, (1.9) 

with fixed momentum p and velocity pj (p). In our context, p will be small for the 
normalized solutions, whereas q can be arbitrary. 

The above scenario is of course the analogue of the center-stable manifolds of j33j 
E] for NLS, where the family of solitons has 8 dimensions, parametrized by the non- 
relativistic momentum in M 3 , the center of mass in M 3 , the gauge in S l and the charge 
in M.+. However, the NLS case of the latter two references is quite different from 
the wave case treated here, since the symmetry group of NLS (generated by Galilei 
transforms, translation, modulation, and scaling) does not mix space and time in the 
way that Lorentz transforms do. In addition, the modulation equations in the NLS 
case are first order ODEs, whereas in the wave case they at least appear on first sight 
to be of the second order. A considerable amount of effort is therefore expended 
here on setting up a suitable framework for the modulational and scattering analysis, 
see Sections [2] and El In particular, we will set things up in such a way that the 
modulation equations are of the first order. 

We can summarize the findings of this paper concerning the trapped solutions in 
Theorem 11.11 as follows: 

Theorem 1.2. All normalized solutions in H e which are trapped by Q as t — >■ oo lie 
on a smooth, codimension- one manifold Ai + (Q) in %, and those trapped by —Q as 
t — » oo lie on Ai + (—Q) = —Ai + (Q). Analogously, those trapped by ±Q ast — >■ — oo 
lie on A4_(±Q). The manifold Ai + (Q) is invariant under the flow of (11. ip in 
forward time, and all solutions starting on A4 + (Q) scatter to the six- dimensional 
manifold generated by Q as t — > oo in the following sense: there exist p^ G M 3 and 
C 1 paths q(t) G M 3 with the property that 

q(t) -> Poo/ (Poo), u(t) - Q{poo, q(t)) scatters to 0, (t -> oo). (1-10) 
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The center-stable manifolds Ai + (Q) and Ai~(Q) intersect transversely in a codi- 
mension two manifold, the center manifold, which characterizes case (9) of Theo- 
rem \l.l\ 

Finally, all normalized solutions u(t) with energy E{u) = E(Q,0), and which are 
trapped by ±Q in positive times, are characterized as follows: either u is constant 
and equal to some translate of ±Q, or for some to G R, Xq £ R 3 , 

u(t,x) = ±W + {t + t ,x + x ) or u(t,x) = ±W.{t + t ,x + x ) (1.11) 

where W± are solutions of (II. ip which approach Q exponentially fast in H, and 
either blowup or scatter to as t — > — oo. An analogous statement holds for the 
negative time direction. 

For more detailed asymptotic behavior of in Theorem ll.2[ see Section [8] below. 
The solutions W± are the analogue of the threshold solutions found by Duyckaerts, 
Merle [HJ [32] for the nonlinear wave and Schrodinger equations with the H l critical 
power. As in the radial case [22], our proof of Theorem 1 1 . 2 1 relies on the gap property 
of the linearized operator 

L + := -A + 1 -3Q 2 , (1.12) 

i.e., spec(L + ) H (0, 1] = without threshold resonance at 1. This property was 
verified in [13] by a numerically assisted Birman-Schwinger method, whereas an 
analytical approach to the study of such spectral problems is developed in [12] . 

The definition of center-stable manifolds given by Bates and Jones [3] does not 
refer to the asymptotic behavior of the solutions originating on them. Rather, these 
manifolds are defined by their local in time invariance, their tangent space at the 
equilibrium, and the fact that they are transverse to the unstable manifold (and they 
do not need to be unique with these properties), see also Hirsch, Pugh, Shub [18]. 
In particular, the equilibrium must be a true stationary solution which is not the 
case in this paper due to the symmetries of (II. ip . This is why [3] only applies to 
the NLKG equation (II. ip in the radial setting. Moreover, since the method of proof 
in [3] relies entirely on energy estimates and does not use any dispersive properties 
of the nonlinear flow, no asymptotic stability properties of the solutions starting on 
the center-stable manifold are obtained. However, using the action J as a Lyapunov 
functional, Bates and Jones do show that their center manifold is orbitally stable in 
both time directions (in particular, solutions starting on the center manifold exist 
for all times). 

In contrast, the point of view of [331 IS] is that of asymptotic stability theory 
of solitons. Since the equations under considerations in these papers, such as the 
cubic NLS equation in three dimensions, are orbitally unstable and in fact small 
perturbations of solitons can lead to finite time blowup, any asymptotic stability 
properties one hopes to recover need to be conditional. Since the ground state 
soliton in these settings is known to have exactly one exponentially unstable mode, 
the condition turns out to be a codimension one manifold. One advantage of this 
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approach is that it is not restricted to the radial setting and allows for moving 
parameters. The latter always correspond to symmetries of the equation and are 
controlled by suitable orthogonality conditions relating to the discrete spectrum of 
the linearized operators. 

Working in an invariant topology such as the energy class, Beceanu [5] was able to 
show that the center-stable manifold obtained by the asymptotic stability approach^ 
has the properties required by [3J. Strictly speaking, due to the underlying sym- 
metries present in the Schrodinger equation (even in the radial setting one has the 
modulation and dilation symmetries, see pOj) or the NLKG equation (II. ip . a direct 
comparison between the manifolds of [3J and those in [5] , [3_U] , or of Theorem 11.21 
above, is inadmissible. However, passing to the quotients by the symmetries does 
allow for a comparison in some sense. For example, W + , W_ and Q form the stable 
manifold associated with Q. 

2. The Complex Formalism 

We consider the nonlinear Klein-Gordon equation 

w-Aw + w = w 3 , w(t,x) : M 1+3 ->■ R, (2.1) 

in the energy space (w(t),w(t)) G H = iJ 1 (lR 3 ) x L 2 (R 3 ). It is convenient to 
formulate the equation in terms of the complexified variable 



u:=Vw-iw, V := y/l - A. (2.2) 

Using the notation 

n 1 :='D~ 1 Reu, u 2 :=Imii (2.3) 
for any complex function u, the conserved energy becomes 
\u 2 \ 2 + |Vwi| 2 + |«i| 2 |w!| 4 - 



E(u) = / ^ - 1^- dx = \\u\\l 2 /2 - || Wl ||l 4 /4, (2.4) 

and the energy (or phase) space for u(t) is the R-Hilbert space defined by 

Sj:= L 2 (R 3 ;C), (f\g) := Re / f{x)g{x)dx (f,g G $}). (2.5) 

We can rewrite ( 12. ip as the system consisting of (I2.3P and 

Ut = iDu — iu\. (2.6) 

In the following, the relation (I2.3P is always assumed for any vector expression, but 
not necessarily (12. 2p . The ground state in this formulation is denoted by 

O := VQ. (2.7) 



3 Which is called Lyapunov-Perron method, see [37] and [2] , in contrast to the Hadamard method 
used by Bates, Jones [3]. 
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We decompose any solution u of (12.61) in the form 

u = (Q + v)(x- c(t)). (2.8) 

Then the equation for v is 

v t = WCv + d k (£l + v)c k -iN( Vl ), (2.9) 

where £ and N are defined by 

£v = v-W~ l Q 2 v l , N(y 1 )=3Qv 2 1 + (y 1 Y = 0(\\v\\l), (2.10) 

respectively. Note that v does not satisfy (12.2]) in general, due to the time-dependent 
translates c(£). However, this will not affect the argument. 

£ is a self-adjoint IR-linear Fredholm operator on fj. The linearized operator iD£ 
has the following generalized eigenfunctions 

zZ>£VO = 0, iV£iVQ = -V£l, iV£g± = ±kg±, (2.11) 

where g± is given in terms of the ground state p of L + = T> 2 — 3Q 2 : 

9± = {2ky l ' 2 Vp t i(k/2)V 2 p, (L +P = -k 2 p, \\p\\ 2 = 1, p > 0). (2.12) 

The natural symplectic form associated with this complex formalism is 

u(u,u) := (iV^ulu). (2.13) 

Indeed, in keeping with the classical Hamiltonian formalism we need to be able to 
write the equation ( 12. 6ft in the form ii = X E {u) where Xe{u) is the Hamiltonian 
vector field associated with E. This means that 

dE{u)(-)=u(;X E {u)) 

or in other words, 

(u - T>~ x u\\-) = u(-,Wu - iu\) 
which yields (l2TT5jl . We have, for a, = 1, 2, 3, 

u(d a £l, idpQ) = S ai/ 3\\dxQ\\l = 8 a ^J{Q), u(g±, g T ) = ±1, 

w(V0,^±) = ui(iVQ,g±) = 0. 
We now perform the following symplectic decomposition of 

« = K9+ + A_y_ + 7, \± := u(v,g T )/u(g±,g T ) = u(v,±g T ), (2.15) 
which implies that 

w(7,0±) = (7i|p) = (7a|p) = 

We can then expand the energy 

2(E(u) - J(Q)) = (L +Vl \vi) + IHI2 - 2C(v) = (£v\v) - 2C{v) 

= -2&A + A_ + (£ 7 | 7 ) - 2C(y) (2.16) 
= 2k(\ 2 -\ 2 ) + {£^)-2C(v), 



(2.14) 
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where Xj and C are defined by 

Ai = ^_^, A ^ = ^T^' C(v) = (Q\vf) + ^i = 0(\\v 1 f Hl ). (2.17) 

3. Parameter choice 

We can reduce the number of coordinates by using the Lorentz invariance. To 
be more specific, let w be a global strong energy solution of ( 12. ip and consider the 
following example of a Lorentz transform of w given by 

w(t, x) \-> w u (t, x) := w(t coshf + x\ sinh u, X\ cosh v + t sinh u, X2, X3) (3.1) 

for i/gR, Then one checks that w v is again a strong energy solution of ( 12. ip and 
E(w„) = E(w) cosh v + P\ (w) sinh u, 

Pi{w u ) = Pi(w) cosh v + E{w) sinh u, P a (w v ) = P a (w), (a = 2,3), 
where P = (Pi, P3) denotes the total momentum 

P(w) = (w t \Vw) = ^uj(u, Vm), (3.3) 

cf. (12 .2p . Note that the Lorentz transform preserves E 2 — \P\ 2 . Every solution with 
finite energy \E\ > \P\ is transformed, by a unique element of the Lorentz group, to 
another solution with zero momentum: 

= P(u) = u(u, Vu) /2 = cj(v, V0) + Vu) /2, (3.4) 

which minimizes the energy among the family of solutions generated by the Lorentz 
transformation group. By the discussion in Section[T]the cases \E\ < \P\ are covered 
by [31] and can be ignored. Since the dynamical properties such as scattering and 
blowup are not changed by the Lorentz transforms, we may restrict our dynamical 
analysis to the invariant subset 

fio ■= {uEfi \ oj(u, Vu) = 0}. (3.5) 

Note that this is well-defined since u gains a derivative. However, we will keep 
the generality of in the static analysis, see Lemma 14.21 below. We define the 
decomposition 

u=s(Q. + v)(x-c), s = ±l, ceR 3 (3.6) 
by the orthogonal projection in H~ 1 (M. 3 ), or the minimization 

||ui -sQ(x - c)\\ L 2 = mm Wm =F Q(x - fe)|| L 2, (3.7) 



which is attained for any u G Sj, and uniquely so if the right-hand side is small 
enough. 

The minimization implies the orthogonality conditions 

= u(v,id j Q) = (v 1 \d j Q) (j = l,...,3). (3.8) 
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Differentiating (13. 8 ft with respect to t, we obtain the parameter evolution in Sjo 

= d t u)(v , id a Q) (a = 1,2, 3) 

= ou(iVCv + c-(£l + v)- iN{v x ), id a Q) (3.9) 

= IIQiQIIiaC - (Vd a Q\ Vl ) ■ c - u(v, Vv)/2, 

where we used (13.41) to rewrite the first term of the second line, while the nonlinear 
term N does not contribute because of its reality. Thus, we obtain 

\c\<\u(v,Vv)\<\\v(t)\& (3.10) 

as long as ||^i||l 2 <^ 1 and the solution u G $) . (13. 8p and the symplectic decompo- 
sition (12.151) imply the orthogonality relations 

= ( ll \d J Q) = ( ll \p) = ( l2 \p). (3.11) 

Since {p, VQ} covers the non-positive eigenfunctions of L + (see [Ml Appendices A 
and E]), 

(Cj\j) = (L +7l | 7l ) + || 72 ||2 ~ || 7l ||^ + || T2 ||| ~ || 7 ||2. (3.12) 
The linearized energy norm \\v\\e > is defined in the subspace (13.81) by 

IMI! = *(A? + A) + mi)/2 = ^(\ 2 + + A 2 „) + <£ 7 | 7 >/2 ~ || v ||£. (3.13) 
In fact, we have 

\\v\\ E ~ inf \\u =F 0(x - b)\\ L 2, ( 374 ) 

since 

\\viWh 1 ~ IK =F Q( x ~ b)\\w < II ± Q( x ~b) - Q(x - c)\\ H i 

~ || ± Q(x - b) -Q(x-c) || i2 

< ||u x =f Q(x - b)\ L 2 + ||ui - Q{x - c)\\ L 2 

< 2\\u l TQ(x-b)\\ L 2 
and so, ||t>i||jf-i < \\ui =)= Q(x — b)\\ H i. Moreover, 

\\v\\ 2 E = E(u)-J(Q) + 2k\l + C(v), C{v) = o{\\vf m ). (3.16) 

We thus obtain the following lemma, which introduces the nonlinear distance func- 
tion from [2H] in the non-radial setting. 



(3.15) 



Lemma 3.1. There exists 5e > and cIq(u) : — > [0, oo) continuous such that 

d Q (u) ~ inf ||tt =f£2(6)||l2, 

±,fc (3.17) 
d Q (u) <5 E d Q (u) = E(u) - J{Q) + 2k\{. 

Moreover, (13 ,7p has a unique solution ($, c) for (Iq(u) < 5e, and decomposing 

u = s(£l + v)(x-c), v = X + g + + A_5-_ + 7 , \± = u(v, ±g T ), (3.18) 
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we have 

d Q {uf ~ |M|! = k -{\\ + A 2 _) + i(£ 7 | 7 ). (3.19) 

In addition, if 

2(E(u) - J(Q)) < d Q {uf < 5% (3.20) 
then d Q (u) ~ |A X | = |A+ + A_|/2. 

4. Hyperbolic and variational estimates 

Next we investigate the hyperbolic structure in Sj . The evolution of A is obtained 
by differentiating f l2~T5j) . Via §ZM) this yields 

A± = u(WCv + c ■ V(£2 + v)- iNfa), ±# T ) ^ 



±[fcA± + (2fc)" 1 / 2 (iV( Vl )|p) - 



and so 



Ai = fcAs - \/ Tr^ilVp) • c, 

V 2 (4.2) 



Both of these equations exhibit the hyperbolic nature of the ODE for (A + ,A_) or 
(Ai,A2). Hence 

d t d Q {u) = 4A;AiAx = 4& 2 AiA 2 + 0(Ai ||u|||) (4.3) 

in the region c?q(m) < Moreover, in ( I3.20p . one has Ai ~ — sdq{u), with s = ±1. 
If in addition d t dg(u) > 0, then 

sAi ~sA+ > sA_ -O(Xf). (4.4) 

Let v g := A + g + + A_g_. Then the projected energy 

E g (u) := J{Q) - fcA+A_ - C{v g ) (4.5) 

solves the equation 
d t E g {u) = -k\ + X. - fcA+A_ - (JV(v)l^i) 

(4.6) 

= (N( Vl ) - N(v gl )\d t v gl ) - u(v, Vd t v g ) ■ c = 0(|| 7 || 2 |M| 2 + \\v\\j), 
in the region (13.201) . and also 

E{u) - E g {u) = (£ 7 | 7 >/2 - C{v) + C(v g ) ~ || 7 || 2 + 0(|| 7 || 2 |M| 2 ). (4.7) 
Therefore, in the region (I3.20p we have 

fLi{o,T) ^ Il7(0)||l + l|7lUfii(o,T)ll^ll(£«ni») t (o,r) + II A II 1^(0,^) - ( 4 -§) 
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Next we compute the leading term in the kinetic functional, defined by 

3, 



m 



dx, K 2 (ip) :- 



\Vtp\' 



Inserting u — 0. + v and using the equation of Q, we obtain 

K (u 1 ) = -2(Q 3 \v 1 ) + O(\\v 1 \\ 2 Hl ), 
K 2 (u 1 ) = -(Q 3 -2Q\v 1 )+0(\\v 1 \\ 2 Hl ). 

| Aip + 71, and L + Q = —2Q 3 , L + p = —k 2 p, we have 



dx. (4.9) 



(4.10) 



Since V\ = 

Mm) 
K 2 ( Ul ] 



k\ T (Q\p)X 1 -(2Q 3 \ ll ) + 0(\\v 1 \\j I1 ) 



(4.11) 



(2 + k 2 /2)\j -(Q\p)X 1 - (2Q + Q 3 | 7 i> + 0(\\ Vl \\ 2 HK 



We can now formulate and prove the following nonradial version of the ejection 
lemma from 1291. 



Lemma 4.1. There exist constants 5x G (0,Se), C* > 1 with the following proper- 
ties: Let u be a local solution of the NLKG equation ( 12. 6 p on [0, T] which belongs to 
£)o and such that 

R := d Q {u(0)) < 5 X , E{u) < J{Q) + R 2 /2, 



(4.12) 
(4.13) 



and so that for some t G (0, T), 

d Q (u(t))>R (0<Vt<t )- 

Then u extends as long as dq{u(t)) < 5x, and meanwhile, 

d Q (u(t)) ~ -sAi(t) ~ s\+(t) ~ e kt R, 

|A_(f)| + || 7 (f)IU» <R+(e kt R) 2 , 
sK s ( Ul (t)) > d Q (u(t)) - C*d Q (u(0)), 

for s = 0, 2 and with either s = 1 or s = — 1. Moreover, dQ(u(t)) is increasing for 
t > R 2 , and d Q {u(t')) > d Q {u(t)) - 0(R 5 ) for0<t< t' < R. 

Proof. We have shown dq{u) ~ — sXi, as long as R < dq{u) < 5e- f !4.13|) implies 
9 t cZg(w) | t= o > 0, and so A + (0) ~ Ai(0). Integrating the equation for A± yields 



(4.14) 



IA- 



3 ±fct A ± (0)| < 



< 



^*-)[|<JV(« 1 )|p>| + ||t;|| i .|c|]( a )i S 
3 fc( *- s) |Ai(s)| 2 rfs, 



(4.15) 



from which by continuity in t we deduce that as long as Re kt <C 1, 

X 1 (t) ~ A+(t) ~ - Sj Re fci , |A±(t) - e ±fct A ± (0)| < (.Re^) 2 . 



(4.16) 
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The estimate on 7 follows from this and ( 14.81) . 

The equation for A 2 together with d t dQ(u)(0) > implies that — sA 2 > R(e kt — 1) — 
0{R 3 ). Hence, for t > R 2 we have -sA 2 > R and d t d Q (u) > 0. For < t < f < R, 
we have d t dq{u) > — R A and so dQ{u{t')) > dqiuit)) — 0(R 5 ). The estimate on K s 
follows from the dominance of Ai, together with (14.111) . □ 

Next, we formulate the important variational lower bound from [29]. Note that 
in this case we work with Sj and not Sj . 

Lemma 4.2 (Variational lower bound). For any 5 > 0, there exist e (5), n , Ki(<5) > 
such that for any u G S) satisfying 

E{u) < J(Q) + s (5f, inf \\uTQ(- + b)\\v > 5, (4.17) 

one has either 

#oM < -ki(5) and K 2 (u 1 ) < -Ki(5), (4.18) 

or 

Ko(ui) > mm(Ki(8) , k \\ui\\ 2 h i) and K 2 (ui) > min(Ki(<5), k \\ Vui||| 2 ). (4.19) 

Proof. The proof is essentially the same as in the radial case [25]. The only major 
difference is that we cannot use the compact imbedding H^ d ^ L^(M 3 ); instead, 
we shall use the concentration compactness method of Lions, [26], [27]. 

We first prove the statement separately for Kq and K 2 , by contradiction. Let 
u n G -f) be a sequence satisfying (I4.17P with e = l/ n but neither (I4.18P nor (I4.19p . 
with either s = or s = 2 fixed. In particular, K s (ui) — > as n — > 00. Since E{u n ) 
is uniformly bounded, we conclude that is uniformly bounded in H 1 . First, 

if || 2 0, then by Sobolev imbedding also ||u™ H4 — > as n — > 00 whence also 
|| Vu™|| 2 —> 0. But then the lower bound in (I4.19P does hold, which is a contradiction. 
Hence, we may assume that ||tt™ || 2 — > Co > as n — > 00. 

We apply the concentration compactness argument to the bounded sequence G 
H l {E?). In the compactness case, there exists a sequence y n G M 3 so that 

«i (• + 2/n) -> ^00 strongly in L 2 n L 4 

as n —> 00. Let 

GoM := Jfa) - = II^H^/4, 

G 2 (^) := Jfo,) - K 2 (<p)/3 = \\Vp\\ 2 L2 /6 + 11^1^/2 
But then Vw™ — ^ Vu^ weakly in L 2 (IR 3 ), and 

G s ( Uoo ) < J(Q), K s ( Uoo ) < 

and 7^ (the latter holds since = would entail that ||Vit"|| 2 — > 0, which 
leads to a contradiction as before). But then Uoo = ±Q(- + b) for some b G M 3 
by the variational characterization of Q (see, e.g. [191 Lemma 2.4]). This means 
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that ||Vu'™|| 2 —> ||Vuoo||2 = HVQIh, whence also it" —> Uoo strongly in H 1 . This 
contradicts (I4.17p . 

In the vanishing case, one has H4 — > 0. But then ||tt"||#i — > 0, which leads to 
a contradiction as before. 

It remains to treat the dichotomy case. Thus, let v n and w n be the separating 
sequences. Since their L 2 norms converge to non-zero values, the positive functionals 
G and G 2 converge to some values in (0, J(Q)), for both sequences. Since Q is the 
minimizer of G s in the region K s < 0, we deduce that K s is positive for v n and w n 
for large n. On the other hand, 

limsup[if s (f n ) + K s (w n )] < 0, 

by the choice of u n , and so K s — > for both v n and w n . Since J{y n ) = G s (v n )+o(l) < 
J(Q), and the same for w n for large n, we deduce, from the lower bound on K s below 
the ground state energy [T21 Lemma 2.12], that v n and w n tend to in H 1 , and so 
in L 4 . Finally, so does which places us back in the vanishing case. 

After obtaining the conclusion separately for s = and s = 2, the remaining proof 
in [29] by connectedness works as well in the nonradial case, where A corresponds 
to V2X1 in this paper. □ 

Just as in the radial case, the above two lemmas allow us to define the sign 
functional. We again restrict to fjo. 

Lemma 4.3 (Sign functional). Let 5s '■= 5x/(2C*) > where 5x > and C* > 1 
are the constants from Lemma \4-l\ Let < 5 < 5s and 

%) := {u E i5 I E(u) < J{Q) + min(d Q (n) 2 /2, e (5) 2 )}, (4.21) 

where Eq(5) is given in Lemma \4-S\ Then there exists a unique continuous function 
<5 : S)( S ) -> {±1} satisfying 



(4.22) 



u€fi(s), d Q (u)<5 E 6(u) = -signAi, 

w G dg(u) > 5 =^ &(u) = sign K (ui) = sign K 2 (u 1 ), 

where we set signO = +1. 

Proof. The proof is the same as in the radial case [29]. □ 



5. One- pass theorem 

Once we have obtained the hyperbolic and the variational estimates, the one-pass 
theorem is proved almost in the same way as in the radial case [29]. The only 
remaining, but minor, differences are 

(1) The center c may be different for an "almost homoclinic" orbit, between the 
departing time and the returning time. 

(2) In the hyperbolic region, the nonlinear distance function (Iq(u) is not strictly 
convex in time. 
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For the issue (2), we have only to ignore small fluctuations around the "bottom" 
of the hyperbolic trajectory. After choosing the small numbers e, R, 5* > as 
in [22], let u(t) be a solution on the maximal interval I C M satisfying for some 

n < r 2 < r 3 e I, 

E(u) < J(Q)+e 2 , max d Q (u(Tj)) < R < R + R 2 < <Iq{u(t 2 )). (5 u 

3=1,3 \ ' 1 

Then there exist Ti G (ri, T2) and T2 G (Y2, r 3 ) such that 

d Q {u(T x )) =R = d Q (u(T 2 )) < d Q (u(t)) < t < T 2 ), 

max d Q (u(t)) > R + R 2 . ^ 

Ti<t<T2 

Let Ai be the totality of minimal points t m G [T 1; T 2 ] of dQ(u(t)) with a minimum 
< <5*. By the ejection Lemma [4.1[ we can extract a finite sequence t\ < t 2 < ■ ■ ■ < 
t n G Ai for some n > 2 such that ti = Ti, t n = T 2 , and for each j = 1, . . . , n — 1, 

tj ^ +i dQ(u(t)) > 5 X , (5.3) 

and for some i,+i/ 3 < t j+2/3 G (tj,t j+1 ), 

'd Q {u(t)) ~ e fc l*-^ld Q (u(t i )) (i,- < t < t i+ i /3 ), 
dg(«(t i+ i/ 3 )) = 5x = d Q (u(t i+ 2/ 3 )), ,g ^ 

dqiuit)) > R* (tj+1/3 < t < t j+2/3 ), 

d Q (u(t)) ~ e fe l*-*^ld Q (M(t i+1 )) (i i+2/3 < i < 

Note that the assumption dq(u(T 2 )) > R + R 2 was used to have at least one time of 
ejection up to 5x between 7\ and T 2 , otherwise u could be just fluctuating around 
dq ~ R on the whole interval [Ti, T 2 ]. 

The modification required by the other issue (1) is also straightforward. We define 
the cut-off function w(t,x) for the virial identity by 

w = X ((x - c(T0)/(t -T l + S)) X ((x - c(T 2 ))/(t -T 2 + S)) (5.5) 

for some S > satisfying |logi?| <C S <C where xO^) £ C^°(]R 3 ) is a fixed 

radial function satisfying x{ x ) = 1 f° r \ x \ — 1 an d X^) = for \x\ > 2. We have 
the localized virial identity 

V w (t) := (wu t \(xV + Vx)u), V w (t) = -K 2 ( Ul (t)) + 0(E ext (t)), (5.6) 

where E ext denotes the exterior energy defined by 

E ext (t) = - I [\u 2 \ 2 + |V«i| 2 + \ Ul \ 2 }dx, 

2 Jx.it) (5.7) 

x g X(t) \x - c(Tx)| > (t-T 1 + S) or \x - c{T 2 )\ > (T 2 - t + S). 

Then as in [29], we have 

E ext (t) < E ext {Tt) + E ext (T 2 ) < R 2 (Ti <t < T 2 ) , (5.8) 
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and so 

V w (t) = -K 2 ( Ul (t)) + 0(R 2 ) (71 < t < T 2 ). (5.9) 

In conclusion, by the same argument as in [29] we arrive at the following no-return 
statement. 

Theorem 5.1 (One-pass theorem). There are constants > such that 2e* < 

R* < $x with the following property: If u G C(I; Sj) is a solution of NLKG ( 12. 6 p on 
an interval I such that for some e G (0, ej, R G (2e, R*] and T\ < r 2 G /, 

E(u) < J(Q) + e 2 , d Q {u{n)) < R < R + R 2 < d Q (u(r 2 )), (5.10) 

then for all t G (r 2 , oo) n I =: V , we have dq(u(t)) > R. 

Moreover, there exist a constant 5* G (0, 6s) (independent of u) such that < 
£o(&*)) and disjoint subintervals I\, I2, ■ ■ ■ C /' with the following property: On each 
I m , there exists t m G I m such that 

d Q {u(t)) ~ e fc l*-*-ld Q («(t m )), min s^ s («i(t)) > - C*d Q (u(t)). (5.11) 

where 5 = &(u(t)) G {±1} is constant, dq{u{t)) is increasing for t — t m 3> i? 2 , 
decreasing for t m — t ^> R 2 , and equals to <5x on 9/ m . For each t G /' \ |J m J m and 
s = 0, 2, we have (£ — 1, t + 1) C I', dq{u{t)) > 5*, and 

/•t+i 

/ minsAT^ui^))^' > i? 2 . (5.12) 
7t-i s=0 ' 2 

6. Dynamics after ejection 

After the ejection, we obtain the blowup in the region & = —1 by the same 
argument (Payne-Sattinger) as in the radial case [29] . 

For the scattering after ejection in the region (3 = +1, we need some small 
modifications. First, we should replace the linear profile decomposition given in [29] 
with the general version in [19]. We can keep the Strichartz norm L^L®, but it is 
more convenient to use the symmetric L\ x norm. It is standard and easy to see 
that these norms are interchangeable for the solution u of f)2.6p . because by Holder's 
inequality 

IKIUfLg < IMIl^lkllJs/^, \W\\l\L% Z \\ u l\\l rL lhl\\ls L e, (6.1) 

and by Strichartz 

||«1 - V l\\ {LT HlnL S t /3 Ll)(0,T) - \\ U 1\\l\LI{Q,T) < ||«l||?|i6 ( o,T) 

< II II II 112 <- 6 - 2 ) 

^ II«i|Il^(o,t)II«iII l 8/3 l? , 

where v denotes the free solution with v(0) = u(0). Hence each of the two norms can 
control the other (with the aid of the uniform energy bound in the region & = +1). 
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Another issue is that the nonlinear profiles for a sequence of solutions in i5o do not 
necessarily have vanishing momentum. Hence, after constructing the single profile 
u° which is not scattering, we need to Lorentz transform it back into Sjq. Here a 
crucial observation is that ||w||z4 L 4( 0jOO ) < oo is preserved by any Lorentz transform. 
More precisely, we have the following result. 

Lemma 6.1. Let u be a finite energy solution of NLKG ( 12.11) on (T, oo), and let u' 
be a Lorentz transform of u. Then there exists T'eR such that v! extends to a finite 
energy solution on (T', oo). Moreover, if \\u'\\ l a ( t> T') < oo then \\u\\ L 4 ^ Too ^ < oo. 

Proof. Let x £ (^(M 3 ) be the cut-off function as before, and let w be the solution 
of NLKG with w(T) = (1 — x( x /S))u(T). If S >> 1 is large enough, then w is 
global with ||w||i4 (ri+3) < oo, by the small data scattering theory. By the finite 
propagation property, u(t,x) = w(t,x) for \x\ > S + \t — T\. Hence u extends to 
this region for t < T. The image of the region 

(T,oo) x M 3 U {(t,x) | \x\ >S+ \t-T\} (6.3) 

by any Lorentz transform contains (X", oo) x R 3 for some T' > 0. Then the transform 
of u solves NLKG for t > T'. After a suitable rotation, we may assume that the 
Lorentz transform is in the form (13. ip . If Uw'Hi* (t>T>) < 00 > then 

00 > }X {t>T') — 1 1 u 1 1 Lf >x (t cosh v-x x sin v>T' and t>T)- (6.4) 

Since the remaining region 

{(t, x) | t > T, t cosh v — x\ sinh v < T', and \x\ < S + \t — T\} (6.5) 

is bounded in space-time, the L\ x norm of u in that region is bounded by the Sobolev 
embedding. Hence |M| l^^t) < oo. □ 

Similarly, we have a local version of the above: 

Lemma 6.2. Let u be a finite energy solution of NLKG (12. ip on a time interval 
I 3 T. Then there is an open neighborhood O of the identity in the Lorentz group, 
such that the transform of u by any g G O extends to a solution in a space-time 
region including a time slab which contains T. 

Proof. Let w be the global solution as given in the proof of the previous lemma. Then 
u extends to/xR 3 U{|x| > S* + |t — T|}, which is mapped to a region containing a 
time slab by any Lorentz transform sufficiently close to the identity. □ 

Therefore the argument in [29] works in the nonradial setting, by using the Lorentz 
transform of u° with momentum in applying the ejection lemma, and going back 
to the original profile u° when using it to approximate the minimizing sequence of 
solutions. Thus we obtain a critical element as well as its precompactness. After 
that the argument is the same as in Kenig-Merle [23], or more verbatim in [19] 
(neither using the radial symmetry). 

Hence, we arrive at the following conclusion. 
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Theorem 6.3. Let < e < e* and let u G C([0,T);S)) be a solution of NLKG (EE} 
on a forward maximal interval [0, T) such that P(u) = 0, 

E(u) < J(Q) + e 2 , d Q (u(t)) > R,, &(u(t)) = ±1 (0 < t < T). (6.6) 

If & = —1, then T < oo. If & = +1 , then T = oo and u scatters to as t — > oo. 
Moreover, there is M G (0, oo) determined only by e such that ||mi||l4 ( t>0 ) < M. 

7. Global dynamics 

Now we investigate some basic topological properties of the scattering and blow-up 
sets. The minimized energy by the Lorentz transform is given by 

E m (<p) = ^J\E{^-\P{ V )\^\gn{E{ V y - P{ V f). (7.1) 
For any e > 0, we define 

Sj<z = {<p e L 2 \ E m (<f) < AQ)+e 2 }, 

%T = W e L 2 I E m {<p) = J(Q) + e 2 }, (7.2) 
fi- £ = {<p <E L 2 \ E m ((p) < J(Q) + e 2 }. 
For cr = ± and any *, we define 

S* = {u(0) G fj* | tt(£) scatters as at — >■ oo}, 

B* = {m(0) G ij* | «(*) blows up in at > 0}, (7.3) 

r; = ^\(5:us:). 

The definition of the "trapped set" T* appears rather awkward. It follows from the 
preceding that any solution in T+ e is forward global, and after the energy-minimizing 
Lorentz transform, it stays close to {±Q(x — c)} cG r3 within distance 2e for large t, 
provided that e < £*. 

It is obvious that S a D B a — and XI = | if G (where X stands for any 
of the three types of sets). The scattering theory implies that S< 6 is open for any 
e > 0. 

To see that Bf £ is open for e < £„, let it be a solution on [0,T # ) with P(u) = 
blowing up at T*, and we proceed as in the radial case. The local Cauchy theory 
implies ||M(t)|| i 2 > \T* — t]^ 1 ^ 2 , and so from the identity 

d 2 \\ Ul (t)\\ 2 L2x = 2[\\u 2 (t)\\% - K ( Ul (t))} > 611^111, + 2\\ Ul \\ 2 m - 8E(u), (7.4) 

we deduce that cy-ui||? 2 = 2(ui\u 2 ) —> oo and K (ui) —> — oo as t — > T* — 0. In 

x 

other words, for any M G (0, oo) there exists T G (0, T*) such that (liiliia) > ^ 
and i^o(wi) < — M for T < £ < T*. Let t> be a solution with u(T) close to u(T) in 
fj, and let w be its Lorentz transform with P(w) = 0. Since P(v) ~ -P(w) = 0, the 
transform is close to the identity, and so by Lemma |6T2| w is a local solution around 
t = T, and moreover w(T) is close to u(T) in Sj. In particular (u;i|u;2)(T) > M/2 
and i^o(wi(T)) < —M/2. By the above identity, (wilu^) is increasing as long as 
Kq(wi) < 0, and the latter can be changed only if cZq(io) < 5*, which is impossible 
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if (wi\w<2) > M/2 3> 1. Hence (wi\w<2) is increasing, as well as ||u>i||z2, and so w 
blows up in t > T by Payne-Sattinger's argument. Then Lemma [6.11 implies that v 
also blows up in t > T, which means that B+ £ is open in $). Hence T+ e is closed. 

We already know that the 9 intersections of (<S+ e , £>+ e , T+ £ ) by (<Sf e , <B< £ , 7I5 £ ) 
are all non-empty containing infinitely many radial solutions, for any e > 0. The 
construction in [2H] works even if we choose a nontrivial dispersive component 
||7(0) IU 2 s. In this fashion, we easily obtain nonradial elements in each set 
(besides those generated by the invariant transforms). 

S< £ , S= £ and Sf are connected for any e > 0, by the same proof as in the radial 
case in [29"] . 

8. Center-stable/unstable manifolds 

In the rest of paper, we prove that the solutions staying forever close to the 
manifold of the ground states scatter to the manifold, using the dispersive estimate 
for the linearized operator. Here we encounter more difference from the radial case, 
due to the parameter freedom. Indeed, the argument is more similar to the NLS 
case in [5U] . 

8.1. Equations. First we need more detailed analysis of the Lorentz invariance. 
With any p G M 3 , we associate the following Lorentz transform 

p =: s9, 9 G S* 2 , s = sinh v > 0, c := ^/l + \p\ 2 = cosh u, 
p 

t := — = 9 tanh v. xq := 9(9 ■ x), Xi_ = x — Xg, (8 1) 

+ \ P \ 2 { ^ 

u p (t, x) := u(ct — sx ■ 9, c(xq — rt) + x±) = u(tc — p ■ x, x + (c — 1)xq — tp). 
In particular, the static solution Q(x) is transformed into the traveling wave solution 
Q{c{x e - rt) + x±) = Q(x - rt + (c-l)(x- rt) e ) (8.2) 

with the velocity r = p/{p), for each fixed p G 1R 3 . Its trajectory is in the 6- 
dimensional manifold in H 1 consisting of 

Q(p, q)(x) := Q(c(x - q) e + (x - q) ± ), (8.3) 

whose vector version is denoted by 

Q{p,q):={V + iT-V)Q{p,q). (8.4) 

The NLKG equation for traveling waves is transformed into the following equation 
of £2{p,q): 

WO. + r ■ VH = tO.1 (8.5) 
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where Q.(p, q)f = [0.(p, q)i] 3 = Q(p,q) 3 (and similarly for other powers). Differenti- 
ating (j8.5p in (p,q), we obtain the linearized equations 

{WC + r ■ V)d q Q = 0, d q Q = -VO, 
(iVC + T-V)d p Q = (d p T)-(d q Q), 
where the M-linear self-adjoint operator £ now depends on the parameters: 

£ = C(p, q) = l- ZV^Qip, q)\V~ x Re . (8.7) 

Note that this agrees with our previous definition (I2.10p for p = q = 0. Decomposing 
the solution u in the form 

u — 0.(p, q) +v(x — q), (8.8) 
we derive the equation for the perturbation v, i.e., 

v t = {iVC p + r(p) ■ V)v + {q- r{p)) ■ V(H P + v)-p- d p Q p - iNpfa), (8.9) 
where Q. p := H(p, 0), C p := £(p, 0) and N p (vi) := 3Q, p ivf + vf . For brevity, let 

A p :=iVC p + Tip)-V, -y:=q- [ r(p(s)) ds, 7r:=(p, 7 ) (8.10) 

Jo 

then the equation can be rewritten in the form 

v t = A p v + ^- V(H P + v)-p- d p Q p - iN P ivx). (8.11) 

The natural orthogonality condition is 

= uiv,d g Q) = u(y,dp£X). (8.12) 

Here, and in what follows, it will be understood that all derivatives of H are to be 
evaluated at q — 0. In particular, (18. 4p implies that 

d p £(0, 0)(x) = zVQ(x), d q Qi0, 0)(x) = -V0(0, 0)(x) = -VPQ(x), (8.13) 

which constitute the root modes for iVC, see (12. lip . By differentiation of (I8.12p . 
we obtain the parameter evolution 

Q = u(v t , d a Q) + uiv , d a d p Q)p 

= -u(v, A p d a O) + w( 7 V(n + v)- pdpO., d a O) (8.14) 

- w(ziV p (vi), <9 a H) + w(v, d a d p 8)p, 

for a = qi,q 2 ,q3,pi,p2,P3- In view of (18. 6p and (I8.12p the first term of (I8.14p 
vanishes, whence 

7[w(9 g n, d a Q) - u(v, d q d a Q)] + p[w(<9p0, <9 a 0) - £j(u, d p d a £l)] 

= -u(iN p (v 1 ) 7 d a Q) 7 

which implies 



7T 



< l|e" N/ V||L (8.16) 
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provided that the 6x6 matrix u(d a £l p , dpQ. p ) is non-degenerate. Certainly it is at 
p — 0, and therefore remains so for small p by continuity. 

Now we look for a bounded global solution v by a contraction argument, for 
a given initial data near O(0, 0) in the center-stable direction. The contraction 
mapping (v,ir) (v,tt) is defined by 

(d t - A p )v - j t ■ V(H P + v) + p t ■ d p Q p = -ziV p (vi), 
d t u(y, d a Qp) = 0, (Va = pi,P2,p 3 , 9i, 92, 9a) 
with the initial constraint 

= u(v(0), d a £l{0, 0)) (Va), tt(0) = tt(0) = 0, 7(0) = 0. (8.18) 

In (I8.17p . (I8.18p . the soliton 0(p, q) is to be evaluated at q = and p(£) which is 
determined by the given path 7r(i). The orthogonality equation is equivalent to 

7 t u;(<9 9 £2p, <9 a p ) + p t u(d p £l p , d a Q p ) 
= 7ta;(u, d q d a £l p ) + p t o;(5;, d p d a 0. p ) - u(iN p (vi), d a O). 
We apply the symplectic decomposition at 7r(0) = 

1 = P + P 1 = P + P + + P_ + P C , P ± £=:A ± s±, P C £ = F, (8.20) 

where P , P± and P c are projections onto the subspaces spanned respectively by 
{WQ, iVQ}, g±, and the rest, cf. ( 12. lip . The equation of v yields 

(d t =F k)\± = u(a ■ Vv + f,±g T ), a := r(p) +7*, 

(ft- A)2-=P c (a- W+/), (8.21) 

/ := -3i(£2j! - Ogjui + it ■ VO p - p t ■ d p 0. p - iN p ( Vl ). 

Note that k,g±, Aq,P c are all time-independent. For any given A_(0), the unique 
bounded solution of A± is given by 



00 



A + (t) = - y e fe ^a;(/( S ),^_)^, 

A_(f) = e- fct A_(0) - f e k ^u(f(s),g + ) ds, 



(8.22) 



while the equation for z is rewritten as 

(d t -A )z = P c (a-VI+g), g:=a-VP v + f. (8.23) 

8.2. A priori bounds. Inverting the 6x6 matrix in (I8.19p . we obtain the estimate 
||7Tt||(L 1 nL° o ) t ^ |Kt||L 1 nL°°INU°°L2 + llAp^OlltLinL^MLi+L-),- (8.24) 
Since P (p(t))v(t) = 0, we have, in any Sobolev space X on M 3 , 

U< lblUHHU< Hx, (8.25) 
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and so 

INU < £ + \\z\\x- (8.26) 

Applying Young's inequality to (I8.22p yields 

||A||(L2 nL oo )t < |A_(0)| + \\f\\(pr\L<*')t{V-+L<x>) a (8.27) 
and applying Proposition 19.81 to (I8.23P (see Section |9~TTI for the definition), 



\\ z \\L™LZnL'jX % |N(0)|| 2 + \\g\\LlL%+LlY, 

X := B~ 5 2 /6 n V u ' 2 L 2 ^ a , Y := D V~^ 2 L 2 ^, 



(8.28) 



for any v, a > satisfying the condition in the proposition, under the assumption 

IMU« < lblU<»NU + ll7*IU«» ^ 5 < 1- (8.29) 
The nonlinear terms are estimated by Holder's inequality as follows: 

||/ - A^(vi)|| (L 2 ni0 o )t(i i +LOO):c < Iblliocll^l^ + ||ff t || L i n L-, 

ii n n n ii~n n n ii~n (8.3U) 

+ ||7ft||LinL°° + lk||| 42 , 

where we neglected higher order terms. Under the assumption 

|A_(0)| + ||z(0)|| 3 <£«l, |N| L i nL oo + \\v\\ stz < BS, (8.31) 

for some big fixed B > 1 and sufficiently small 5 > 0, we therefore obtain 

IMkaL- < (55) 2 < 5, 
l|u|Uf^ni?x J$ 5 + {B5f < B5, 



(8.32) 



by a bootstrap argument. Moreover, those global bounds imply that, as t — > oo, 7r 
converges, A± — >■ 0, and z(t, a; — b(t)) scatters, where 

b(t) := / a{s)ds. (8.33) 

Hence Pov(t) — >■ strongly, and so a; — 6(t)) also scatters, where b = a = 
r{p) + 7t — j- r(p (X) ) converges, although we cannot generally approximate it by 
a (Lorentz transform of a) free solution, since b(t) — tr(poo) does not necessarily 
converge. 
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8.3. Difference estimate. The above argument does not apply to the difference of 
two solutions, due to the term a ■ Vz in the H equation, which causes two problems: 
possible growth in time and derivative loss for the difference. Hence we employ a 
rather weak norm for the difference, with an exponential weight in t and values in 
H~ l . This is still sufficient, because the main issue in the difference estimate is to 
control A+(0) in the contraction argument. To be more precise, fix p > such that 

< 5 1/A < p < min(l, k), (8.34) 

and define the function space G on (0, oo) by the norm 

||/|| G :=supe-^|/(t)|. (8 . 35) 

t>o v ' 

We estimate the difference of the mapping of two given (v J , vr J ) (j = 0, 1), denoted 
by <X = X 1 — X°, in the following norm 

\\<^t\\a + ||<A|| G + IH^IIg^ 1 - ( 8 - 36 ) 
The difference of (18.191) yields 

ll^tllo ^ ^IHh^lU- 1 + |<H + \<^t\ + \<&t\ + \\<v\\h-i-\\g 

/g 

^ dW^WctH- 1 + &P~ x \\<tot\\G + 5p~ 2 \\<7T t \\ G + Sp^WovWa^-i. 
The difference of (I8.22p gives us via Young's inequality 
~ 3 

\\<M\g < |<A_(0)| + -^[\\<v\\ GtH -i + ll^tllc + p _1 ||<7r t || G + \\<v\\ GtH -i}. (8.38) 

For <z, we need a change of variables to avoid a derivative loss due to the transport 
term. Let b(t) be as in (18. 33 p . and define 

r b <p(x):=<p{x-b{t)), ((t) = T b z(t), ({t)=nz{t), T b A = T b Ar b -\ (8.39) 
where A is any operator on L 2 (IR 3 ). Then (I8.23P is transformed to 

(d t - Ao)( = [3i(T b Q 2 - Q 2 ) - T b P d a ■ V]C + r b P c g. (8.40) 
Now we employ the linearized energy of regularity level H^ 1 : 

IHH-i := {L{iVL)- x P c i ] \{iDLY 1 P c i ] ) ~ \\{%DC]- x P c r]\l - \\Pj)f H -u 

fitlhlH-i = 2(£(iV£)- 1 P c (d t - WC^KWC^P.T)) (8.41) 
<\\{d t -W£)P c r}\\ H -i\\r)\\ E -i. 

From the difference of (I8.40p we infer that 

d t \\<\\ E -i < 6[\<b\ + |<tt| + |<5r t | + Klk-i + \\M\h-A, ( 8 - 42 ) 



and integrating it 

\G t H, 



<\\ Gt H-i £ ll^(0)|| 2 + S[p 2 ||<vr t || G + p 3 \\<7T t \\ G + p 1 ||<r;|| G . fl -i]. (8-43) 
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The difference before the translation is bounded by 

\\<n H ^ - IKIIh- < \\(no-T bl )z°\\ H ^ < \<b\\\z°\\ Ll <6\<b\, (8.44) 

and so 

ll^lUi^- 1 ^ \HWctH- 1 + ^[p~ 2 IKJg + P" 3 ||<^||g]- (8-45) 

Thus we obtain 

IH^Hg + \\<w\\ GtH -i < |<A-(0)| + ||<z(0)|| 2 + 5p- 3 \\<7r t \\ G + *p _1 ||<T;|| Gtfls -i. (8.46) 

In conclusion, there is a unique fixed point (v, it) in the set (|8.31j) . where we have 

b = a = r(p) +7 = q. (8.47) 

It depends continuously on the initial data (A_(0), z(0)) G 1R x L\ in the above 
topology. The same estimate holds for the derivatives with respect to the data, which 
implies that v(0) of the fixed point is smoothly parametrized by (A_(0), z(0)) G 
M x P C (L 2 ). We have already seen in Section 18.21 that, as t oo, X(t) — > 0, 
7r(t) 37TOO, b[t) = a(t) = q(t) r(p 00 ), and 

v(t, x - bit) + 6(0) - q(0)) = u{t, x) - Q(p(t), q(t))(x) (8.48) 

scatters. This means that, at least in a weaker or localized topology, the solution u 
converges to the family of ground states. Moreover, we have 

E(u) = E(£1 P J + \\ Voo \\l/2 = J(Q)( Poo ) + |M|^/2, 
P(u) = P(0 p J + P( Voo ) = J(Q) Poo + P( Uoo ), 

where := lim^oo e~ wt [u(t) - Q(p(t), q(t))]. 

Note that we started with a fixed parameter at t = 0, normalized to (0, 0) by a 
Lorentz transform, and ended up with some non-zero but small parameter tt^ at 
t = oo. It is in general difficult to reverse this process starting from t = oo because 
of the growth in the modulation parameter, unless working with localized dispersive 
data, which would yield better asymptotic control on the parameter. 

8.4. Uniqueness. Next we prove that any solution which stays forever close to the 
family of ground states is necessarily a Lorentz transform of one of those constructed 
above. Thus, let u be a solution of NLKG satisfying 

sup inf \\u{t) - £}(x - q)\\ 2 < 8. / g 5q n 

First notice that this property is preserved by Lorentz transforms which are 0(5) 
close to the identity, provided 5 > is small enough. This is because the above 
condition implies that for some R{5) > and at each t > there is a ball of radius 
R in IR 3 such that the linear energy of u is at most 0(5 2 ) in the exterior of the ball. 
Then by the same argument as in Lemma 16. 1[ the solution extends at least to the 
exterior of the light cones emanating from this ball, with a uniform energy bound 
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of 0(5 2 ), while the interior energy can be controlled by the linear energy inequality 
for time 0(5), whence the claim. 

Next, the implicit function theorem implies that there is a unique (p, q) £ M 6 such 
that for u p (0, x — q) = + f (0, x), 

\p\ + \H0)h £ 5 > = w(«(0), <9 Q £2). (8.51) 

Hence replacing it by u p (t, x — q), we may assume, in addition to (18.50 p . that 

u(0)=jQ + u(0), u(v(0),d a Q) = 0. (8.52) 

Let 7r(t) = (p(t),7(t)) and u(t) = Q(p(t), q(t)) + v(t,x - q(t)) with q(t) = 7(t) + 
Jo t (p( s )) ds where ir is evolved by (I8.14p . which is locally uniquely solvable as long 
as |vr(t)| + ||t'(t)||2 «C 1, preserving the orthogonality u(v(t),d a Q.(p(t),q(t))) = 0. 
Meanwhile, (I8.50P implies, via the implicit function theorem, that at each t there is 
a unique (p(t),q(t)) £ M 6 of size 0(5) such that the orthogonality property holds. 
This gives an a priori bound of size 0(6) on 7r(t) solving (I8.14j) . so it extends to all 
t > 0, satisfying 

IM|l°°(o,oo) + \\v ||i-(o,oo;Ll) < 5, = d a Q(p(t), q(t))) . (8.53) 

To establish the global Strichartz bound, consider the local version for any T > 

Nt '■= \Wt\\L\{Q,T)nLf{T,oo) + IM|L2(0,T;X)nLj?°ffi- (8.54) 

Since (n, v) is bounded, we have (I8.19p . (I8.22p and (I8.23P with (7?, A, z, v) = (tt, A, z, v). 
Then by the same argument as in Section [8721 together with (]8.53p . we obtain 

No<S, N T <8 + 5N T + N*, (8.55) 

uniformly in T > 0. Hence by the dominated convergence theorem, we conclude that 
^oo ^ 5. Now that (n,v) is a solution belonging to the set (I8.3ip . the uniqueness 
follows from the contraction principle via the difference estimates of Section 18.31 

8.5. Threshold solutions. Finally, we prove the Duyckaerts-Merle type classifica- 
tion of the solutions with the threshold energy E(u) 2 — |P(n)| 2 = J(Q) 2 , that is, 
there are only three solutions modulo the symmetries. 

Let u be a solution with E(u) 2 — |P(w)| 2 = J(Q) 2 , satisfying (I8.50p . Then after 
a Lorentz transform, it is given by the above construction, and (I8.49P implies that 
Woo = if E(u) 2 - \P(u)\ 2 = J(Q) 2 . Hence z(t) = P c v(t) -> strongly in L 2 as 
t — > 00, which allows us to apply the Strichartz estimate, Proposition 19. 8} to z 
starting from t = 00. Consider a weight function 

u T (t) = min(e At{t - T) ,l), (8.56) 

for any fixed fi £ (k/2, k) and T — > 00. Then by the same argument as in Section IS72| 
we obtain 

II ^T^t || L 1 nL°° + \ Wv\\ LrH i nL 2 X < e-» T 5, (8.57) 
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where the right-hand side comes from \\i*re M *A_(0)||x,i n L° o - Taking T — > oo, we 
deduce that 

||e At *7r t || L i nL ^ + ||e^|| LrH i nL 2 X < 5, (8.58) 

which allows us to bound the difference of any two such solutions from t = oo. Let 
(V, 7r J ) (j = 0, 1) be two solutions enjoying the above exponential decay. Then by 
the same argument as in Section 18.31 we obtain from the difference of (18.1 9ft 

IKJl- < <5||e -/xt [|M|tf-i + |<7r t | + |<tt 

^ ^IMIl-h- 1 + IH^IM' 

where S comes from the norm (I8.58P for both (t> J , 7r J ). Similarly, we obtain from the 
difference of (I 



59) 



||<A|Uf> ^ l«A_(0)| + *[||<7r t |Uoo + \\<v\\ L?OH -i}. (8.60) 
From the difference of (I8.40p . we have the energy inequality 

-dtKIU-x < Se'^bl + |<tt| + |<rr t | + \\<\\ H -A- (8.61) 
Since ||C J '(^)||2 as t -> oo, we can integrate the above from t = oo, which yields 

IKIIl-h- 1 £ ^[||<tt*IU?° + \\<4l~hA (8.62) 
The difference before the transport is estimated similarly 

- IKIIl^J 1 £ l|e~^||e^2 || Li ||L- < *||<7rt|U r . (8.63) 



Thus we obtain 

||<7rt||£- + \\M\l T h-^ ^ l<A-(0)|, (8.64) 

which implies that any solution satisfying ( I8.53P and E(u) 2 — |P(u)| 2 = J(Q) 2 is 
uniquely determined by A_(0) G K. of 0(5). Moreover, these conditions are invariant 
under the forward time shift, and A_(t) — > as t — > oo. 

Therefore, all such solutions are classified into three distinct cases: A_(t) > 
decreasing to 0, A_(i) = 0, and A_(t) < increasing to 0. Moreover, solutions 
in each class are translations of the others in the same class. The solution with 
A_(t) = is the exact ground state u = Q, and applying the one-pass theorem 
backward in time, we deduce that the solutions with A_(t) > blow up in finite 
time t < 0, while the solutions with A_(£) < scatter as t — > — oo. 

9. Linearized dispersive estimate 

This section is devoted to establishing the dispersive bounds needed in the con- 
struction of the center-stable manifold. As in Beceanu's theorem j5] on energy class 
(or better) center-stable manifolds for the cubic Schrodinger equation in M 3 for gen- 
eral data, we are lead here to a linear equation involving a small perturbation in 
the form ia(t) ■ V, which reflects the translation invariance. A general approach 
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covering this case was developed by Beceanu in [I] for the Schrodinger case in R 3 , 
and we follow his method in a wide sense. However, the technical details here are 
quite different from those of [I] . For a related, but much simpler, statement see also 
Appendix B of 



9.1. Linear matrix operator. In order to use the Kato theory, we need to make 
our operator C-linear. Define a matrix operator £ on L 2 (R 3 ; C 2 ) by 

such that i£ is an extension of iDC in the sense that 

iV&p = $ ^ i£ (f) = (?) . (9.2) 



<pj \tp 

In the same way, the scalar operators % and T>C are extended to ia^ and a^£. Here 
03 = ( ^ ] as usual. The inner product and the symplectic form u> are extended 



-\ / 
in the following fashion: 

I r 2 

(u,v) := - / } Uj(x)vj(x)dx, fl(u, v) := (a 3 T> 1 u, v). (9.3) 

For any vector u G L 2 (R 3 ; C 2 ), we define scalar functions it ± by 

u + := V~ l {ui + u 2 ), u~ := Ui — u 2 - (9.4) 

Then we have 

£^ = ^ L + if + = ip~, tp~=ip + , (9.5) 

which enables us to determine the spectrum of £ from that of L + , cf. [T3] : 

a(i£) = i(-oo, -1] U 00) U {-k, k} U {0}, (9.6) 

and the complete set of (generalized) eigenspaces of the nonself adjoint operator i£ 
is spanned by 

i£gt = ±kgt, *£V<3 = 0, i£R$ = -V<3, (9.7) 

where 



These functions enjoy the following orthogonality properties: 

o = n(gt,R$) = n(gt,v$), 



(9. 



(9.9) 
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Hence the projection P = P c onto the continuous spectrum is given by 
P c = l-P d , P d = P + P+ + P-, P ± u:=n(u,±g^)gt, 

p u -.= JiQ^Mu, R~i})v~$ - n(u, Y~d)i}~d . (9 ' 10) 

The symplectic orthogonality is rewritten in the scalar form 

= tt(u, g}.) = n{u, gt) ^ = (p, u~) = (p, u+>, 

o = n( u , vq) = n( u , Sq) <=^ o = (vg, «-) = (vg, « + ), 

where (•, •) denotes the standard inner product on L 2 (M 3 ;C). 

9.2. Dispersive estimates for the scalar operator. First we recall the weighted 
L 2 estimates for the resolvent and the propagator of the scalar Klein-Gordon (V) 
with a potential. All the estimates in this subsection are classical, but we provide 
proofs for the reader's convenience. 

Let Bp q (M. n ) denote the inhomogeneous Besov space based on L p (R n ) for any 
n > 1, s G R and p,q G [1, oo]. For brevity, we use the standard notation 

H S :=B^ C s :=^ i0O . (9.12) 

The homogeneous versions are denoted by B^ q , H s and C s , respectively. For s G 
(0, 1), we have the equivalent semi-norms by the difference (cf. [5]) 

\\<P\\b» - II SUp \\ip(x) - <f(x - y)\\ L p \\ L i(da/a)- 

The weighted L 2 space L 2,s (R n ) is defined by the norm 

IMU 2 .* = ll(«>VIU a (K») (9-14) 

for any s G R. Hence L 2,s is the Fourier image of H s . 

We start with the resolvent estimate on the free part £o = which is a version 
of the classical limiting absorption principle [T]. 

Lemma 9.1. In any dimension d > 3, 

sup || (£q - z)"V||l2.-i < IMU 2 ' 1 - (g 15) 

As the proof below shows, we can even place the homogeneous weight \x\ s near 
x = for s < 1/2, but the above weaker estimate is sufficient for our purposes. 

Proof. Since a^D is diagonal, it suffices to consider scalar operators. Let x = r6 be 
the polar coordinate, \ £ Co°0^) radial, x{ x ) = 1 on 1/2 < \x\ < 2, and x( x ) = 
on |x| < 1/3 or |x| > 3. Define operators Sj for j G Z by 

Srfrti) = X(kl)^(2^). (9.16) 

After such a cut-off, we may regard the Besov norms as being defined on R. 
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First we derive (19. 15ft from the following three estimates: for any p > 4, 

SUp llV^Vll^ys < \\<Pr\\l£, SUp \\Sf d / 2 - l ' 2 V \\ Ll{K/l)r < Mtf/*, (9. 17 ) 



and for any p < oo, 



sup 



f(r) 
(r) — z 



dr\ < 



2 

^[||r-VlUi +sup||r 1 -^ i /ll ( ^ ) 1- (9-18 
tit ^' eZ p 



Indeed, we can apply the latter estimate to 
|<(2>-z)-V 



</S 



(r) — 2 



(9.19) 



by setting /(r) := j s (p(r6)ip(r6)r d 1 d6. The L 1 part is bounded by Cauchy-Schwarz 
and Hardy's inequalities: 

\\r- k fhi < \\\t\-*0\U\\\Z\-t$h £ W\\h^U\\h^ < IMIl^IMLw ( 9 -2o) 

since d > 3. Since B x Jl is an algebra, the Besov part is bounded by 

V ||5 h r( rf - fc )/Vll L2 ^/ P ||5,r( rf - fc )/ 2 ^|| L ^ 1/p , 

^ — ' ^e^p.i « p,i 

li-/v|<i 



Therefore, using; fl9~T7|) and C 1 / 2 C B^f , we obtain 

which implies ( 19.15jl . 

To prove (I9.17p . one verifies that for any r 2 > r\ > 

»r 2 



lb(r«| 



¥>r(r0) drl < (1 - ri/^) 1 / 2 ^ 1 -^ 2 !!^^)!!^, 



(9.21) 
(9.22) 

(9.23) 

(9.24) 
(9.25) 
(9.26) 
(9.27) 

Then the real interpolation between (I9.26P and (I9.27P gives us for < /3 < 1, 

Sy/2-(i-^ . -> (9.28) 



by Cauchy-Schwarz. Square integrating in 9, we obtain the first estimate 

II V /2 ~Vll iic y 2 < IMUi < |MI#- 
The complex interpolation between (I9.24p and the trivial estimate 

Ws^Ml^ < \\<p\\h, 

yields a uniform bound for j e Z and a G [0, 1]: 

Combining (I9.25P with the Sobolev embedding, we have another bound 

Sjr d / 2 : L 2 X — > LqL 2 C Lj(By a J^ ) r . 
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Choosing < a < 1/2 and (3(1 — a) = 1/2, we obtain the second estimate of (19. 17ft . 
To prove ( I9.18p . let z = r + is with r, e E M and e ^ 0. If r < 1, then on r > 0, 



|«r> - < «r> - I)" 1 = <r>/r 2 



(9.29) 



which gives rise to the L 1 part of (I9.18p . If r > 0, then there is a unique s > such 
that (s) = r and the integral has a singularity at r — s. But it is easy to see 



r < s/2, 3s/2 < r 
s/2 < r < 3s/2 



\((r)-z) 



- 1 ! < 



'r)/r 2 



=► |((r> - z)- 1 - (s(s)-\r -s) + is)' 1 ] < (s)/s 2 . 

Hence it suffices to bound 

f(s + a)dr 



s/2 



f(s + a) , f(s-a) 



a + is 



\a\<s/2 cr + ie 
where the integrand on the right equals 

f(s + a)-f(s-a) , -2%e 



-a + is 



dt. 



(9.30) 



(9.31) 



cr + is 



a- 



f(s-a) 



(9.32) 



To bound the first term, we employ an integral identity for any function F(t) 



F(t)dt 



2da 



J a j2 



(cr + ief 



IE 



t + ie (t + i£)(2t + i£) 



which follows simply from Fubini. Hence (I9.18P equals to 



s/2 



[f(s + t) - f(s + t- 3a/2)]dt- 



2da 



a/2 



a + is)'- 



F(t) dt, 



+ R, 



where the remainder R is bounded by 



\R\ < 



s/2 



<J 2 + E' 



\\f{s + a)\ + \f{s 



da < 



oo ^ 



(9.33) 



(9.34) 



(9.35) 



while the leading term in (I9.34h is bounded by using Holder in t and then by the 
difference norm of the Besov space 



< 



'/ 2 da 
a- 1 / p \\f(t)-f(t-3a/2)\\ L ,-< 



a 



This finishes the proof of (I9.18p . 

Next, we transfer the estimate (I9.15P to the perturbed operator £. 
Corollary 9.2. The operator £ from (19. ip satisfies 

sup ||(£ - z)" 1 P c (p\\ L 2,-i < \\if\\ L 2,i 



(9.36) 
□ 



(9.37) 



30 



K. NAKANISHI AND W. SCHLAG 



Proof. Since Q 2 T>~ 1 is compact L 2,_1 — > L 2,1 and P := P c = 1 — Pd is bounded 
on L 2 ' -1 and L 2 ' 1 , the compact perturbation argument implies that the resolvent 
estimate (19.371) can be violated only if £ exhibits a resonance. To be more specific, 
let R*(z) = (£* - z)~\ £ = £ + K and X = L 2 ~ l . Suppose (19T3TD fails, then there 
are z n G cr(£) c and G L 2 such that 

||^n||x* 0, ^ n := R(z n )P c ip n , \\ip n \\ x = 1, ip n ^ ip vnw-X. (9.38) 

Applying the resolvent expansion, we have 

ip n = R (z n )P<p n + R Q {z n )K,ip n . (9.39) 

Since P : X* — > X* bounded and Ro(z) : X* — > X uniformly bounded, the first 
term on the right is vanishing. For the second term, one has JCPip n — > JCPip in X*, 
because K. : X — > X* compact. If \z n \ — > oo, then we infer from 

R (z n ) = -z?{£o -z n - 2 )R (z n ) = z-\-l + £ Ro{z n )), (9.40) 

that Ro(z n )f — > in X for any / G S C X* dense. This with the uniform bounded- 
ness implies that Ro(z n ) — > strongly as operator X* — > X. Then ijj n — > strongly 
in X, contradicting (I9.38p . Hence z n is bounded, and so extracting a subsequence, 
we may assume that z n — > 3z G C. Let / n := ip n — R (z n )P(p n . Then f n —±ip weakly 
in X C 5', and so in S' we have 

(£ " Zn)fn = /C^n ~> (£ " = (9.41) 

Hence £ip = zip in the distribution sense, whereas ip = Pip G X, since P : X — >■ X 
continuous. If ip — 0, then /C^ n — >■ strongly in X* and so does ip n — > in X, 
contradicting (19.381) . Finally, using ( 19. 5p . one obtains 

L + ip + = zip" , ip~ = zip + 

which implies that ip + ^ and L + ip + = z 2 ip + . If ip + G L 2 then ip G L 2 is also an 
eigenfunction of £, which is a contradiction. Hence ^ + G X \ L 2 . It now follows 
from the theory of scalar Schrodinger operators that the only possibility would be 
that z 2 = 1 is a threshold resonance of L + . However, in our case this is known not 
to occur, see ITS], [T2]. □ 



By the Kato theory [21], (I9.37P implies the weighted L 2 estimates on the evolution. 
Lemma 9.3. For any e > 0, we have 

(9.42) 



2r2,-l ~ 



Proof. Define an operator T by Tip := (x)~ 1 e ttij (p. Then 

/■oo 

T7= / e- fa£ (x)-V(s)ds 



< II ^-2,1 
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Our goal is to show that T \ L? x ^r L? tx which is the same as TT* : L\ x — > Lj x . Now 

POO 

(TT* f)(t) = / (xy^-^ix)- 1 f(s) ds, t > 
Jo 

where we first (and without loss of generality) consider positive times. Denoting the 
Fourier transform in time by F, one checks that 

(TTT*f){r) = [ (x)~ 1 R(r-iO)R{a + iO)(x)~ 1 (Tf)(a)da. 

By the resolvent identity, the right-hand side is the same as 

x (x) (Ff)(a) da. 

T — O — i(j 

By Plancherel, the boundedness of the Hilbert transform on L 2 , and Corollary 19.21 
we conclude that TT* is bounded on Lf x , which proves the first inequality in (I9.42p . 
Similarly, the retarded estimate on t > follows by letting e — > +0 in 



0<s<t 



/CO 
e isT R(r -ie)Pf{s)ds 



\LIX 



r>oo 

<||/ e^f{s)ds\\ L2x ,^\\f\\ LhQX ». 



o 



The estimate on t < is obtained by reversing the sign of ie. □ 

9.3. Time-dependent, spatially-constant perturbation. We now turn to the 
main goal of this section, i.e., proving the dispersive estimates required in the proof 
of the center-stable manifold theorem. To be more specific, we consider a small 
time-dependent perturbation in the form 

7 t = iZ>£ 7 + P c [a(f)-V7 + /], 7 (0) = P 7 (0), \\a\\ Lr < 1, (9.43) 

and follow Beceanu's approach [3] in order to establish Strichartz estimates. By the 
time symmetry, we may restrict it to t > 0. Putting u = (7,7) and F = (/, /), it is 
rewritten in the matrix formulation 

Ut = i£u + P c [a(t)-Wu + F], u(0) = P c u(0). (9.44) 

Let A(i)(p := a(t) ■ Vy. The equation can be extended to include the component 

z t = i£P c z — P d z + A(t)P c z + F, z(0) = u(0), (9.45) 

then u := P c 2 satisfies the previous equation. The equation can be rewritten in the 
form 

z t = [i£ + V + A(t)]z + F, 

~ (9.4b) 
V := -i£ Q P d -P d + iJC, F(t) := -A(t)P d z(t) + F(t) 
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where £ = £o + JC- Let U(t, s) be the evolution operator defined by the equation 
u t = A(t)u, namely 

U(t,s)<p = (p(x + b(t,s)), b{t,s)= [ a{r)dr. (9.47) 



Note that £ and A(t) commute, so do their evolution operators. Regarding F as 
an input, we consider the Duhamel formula 

z = z + [ e it - s)i2o U{t,s)[Vz{s) + F{s)}ds, z := e iZot U{t,0)z{0). (9.48) 
Jo 

Now define space-time operators T,T ,Ti on (t, x) G [0, oo) x R 3 by 

Tg= f e^- s ^ Q U(t,s)Vg(s)d Sj 
Jo 

T g= [ e {t - s)i2o Vg(s)ds J T\g — I e l*-*A»w ) Vg (s) ds 



Then we have 



(9.49) 



(1-T)z = z + / e {t - s)iZo U(t,s)F(s)ds, 

t Jo (9.50) 

(T -T)g= [ e^-^(l-U(t,s))Vg(s)ds 
Jo 

as well as the following properties. 

Lemma 9.4. T ,Ti are bounded on L\L 2 ^~ U for any a > 1, where it holds 

(1-T )(1 + T 1 ) = (1 + T 1 )(1-T ) = 1 
Proof. Let z£ = i£ + V, then 

T TJ = [[ e^*^Veto- tl Wvf(t 1 )dt 1 dt 
J Jo<t 1 <t <t 

= 11 * e^-'^Ve^VfitJdsdt! (9.51) 
Jo Jo 

Jo 

where in the third step we used the Duhamel formula 

e *& ( p = e itto (p+ f e {t - s)iZo Ve isl <fds. (9.52) 
Jo 
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Since i£ = i£ — V, we obtain TiT = — T + 7\ in the same way. The weighted L 2 
estimate f)9.42p for £ implies that if a > 1 then 

\\Tif\\ LfL i,-« < || j\e^- s) P c Vf{s) + e-^e-^°P d Vf(s)} ds\\ L , L ^ 

< II PcVf\\ LlL2A + \\P d Vf\\^,-i < ii/ii^.-, 

since : L 2, ~ <J — > L 2,1 bounded. The bound for To is simpler, as it uses the free 
analogue of (19151) . □ 

Next, we extend the weighted L 2 estimate to T with some small loss of regularity. 
For that purpose, we introduce the bilinear form in dyadic frequency 

poo /•max(t— M.O) 

If (/,<?):= / / (e^ w U(t,s)A jW f(s)\wg(t))dsdt, (9.54) 

Jo Jo 

for j > 0, M > 1, w(x) := f(t,x),g(t,x) E L\ x and 1 = £Jo A ; is a 

Littlewood-Paley decomposition defined by Aju = J r ~ 1 ip(2~^)u(^) for all j > 1, 
with some radial non-negative <^ G (^^(M 3 ) supported on 1/2 < |£| < 2. 
To sum Ij over j > 0, we use the almost orthogonality 

Ij'd-U) Ij'iQjf-Q^l)- Qif= E ^kf + w-^Mf, (9.55) 

|fc-il<i 

where the commutator term is small in the sense 

\\w~ 1 [A k ,w]f\\ 2 <2' k \\f\\ 2 . (9.56) 
For each if , we now have the following estimate. 
Lemma 9.5. Let a > 3/2 and M > 1. T/ien /or any j > 0, 

|/f <J,g)\ < C(a)2*M^f\\ L *Jg\\ Llx , (9.57) 
provided that ||a||L^ <C 1. 

Proof. Define space-time functions Kj and Ri by 

i^(t, x) = 2 3j K j (2 j t, 2 j x) = e iW 2 3j fi(2 j x). (9.58) 
Then we can rewrite If as follows 



If (f,g) = Jf (f,g) + Jf(gJ) (9.59) 

Jf(f,g)= / ,, . Kj(t - s,x - y + b)w(y)f(s,y)w(x)g(t,x)dsdtdxdy 
J J \y\>m, 

s>0, t>s+M, 

2^K\2^u,2h)w(y-b)f(s,y-b) 

\y-b\>\y+z\, 
s>0, u>M, 

x w(y + z)g(s + u,y + z) dsdudydz, 
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where b = b(t, s) = a(a) da. Now we use the pointwise decay of for t > 1, 



{t-^t-N)- 1 / 2 (|t-N|^2^t), (9 - 60) 

where N G N is arbitrary. Note that the first line is similar to the wave equation, 
while the second term is the stationary phase part due to the mass term. The proof 
of (19.601) will be given below. In particular, (I9.60P implies the uniform estimate 

\K j (t,x)\ <2 j rl (9.61) 

whence 

\jf(L9)\<2% [ u-Hw(x)f(s,x)\\ Ll \\w(y)g(s + u,y)\\ Ll dsdu 

Ju>M (9.62) 



<2^M^n Ll Jg\\ Llx 

This bound is not too useful for our purposes since we require more decay in M. 
Using (I9.60p and arguing as in (19. 62ft . the integral (19. 59ft onu^ \z\ is bounded by 



2 3j (2 J 'u)~ 3 ~ N \w(y - b)f(s,y- b)w(y + z)g(s + u,y + z)\dzdyduds 
u>m (9.63) 

<2-^M- 2 - w ||/|| L? J|^|| iL . 
If u ~ \z\, then 

\b(t, s)\ < \t - s\ = u ~ \z\ < \y - b\ + \y + z\ + \b\ < \y - b\ (9.64) 

since \y + z\ < \y — b\. The integral (19. 59ft over this region is therefore bounded by, 
with b = b(s, u), 

2^ \f(s,y-b)g{s + u,y + z)\ 



\u\i\y- b \ 21u (y-b) a (y+zy 



dydzduds 



< / 2 2 >u- 1 \\(y-br°(y + zr°\\ Lf L ,\\f(s,-)\\ L 2\\g(s + u r )\\ Li dsdu 

•/ ii -A I 

< [ 2^u^\\f(s,-)\\ Ll \\g(s + u r )\\ Ll dsdu 

Ju>M 



'u>M 

< 2 2j M^~ a 



\LlJ\9\\L» 



and the lemma is proved. □ 

Proof of ( 19.60p . From the explicit form of 052, one has (up to some nonzero multi- 
plicative constants) 

K j (t, r) = [ e it{p) i +icrp ij(p) dpdc = r" 1 / j t{p) i sin(rp)^(p) dp, (9.65) 
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where (p) . := V ' P 1 + 2 ~ 2j and ip e C£°(l/2,2). If r < t/2, then we use the first 
form. Let $ c (p) := t(p)j + cr P; then 



$ , c = t 7 ^- + cr = 0(t) > & e =t^ = 0(12-*)>\Q£ + Q\, (9.66) 

for all k > 0. Integration by parts in TV times gives 

e'*^(p) dp = / e^z^'J- 1 ]^) ^, (9-67) 



where the right-hand side is bounded by i — ^H^llvv^' 00 - ^ remains to estimate in the 
region r > t/2 and \t — r\ > 1, for which we use the second form in (I9.65p . It can 
have stationary phase only at p = po solving 

2~ 2 H 

= ^(po) = t-r- (9.68) 

{Po) j {(po) j + Po) 

which has at most one solution in supp?/' because <E>" > 0. If \t — r\ ^ 2~ 2 H, then it 
has no solution in supp ip and 

1$'^! ~ \t - r\ + 2~ 2 H > \$ { l? k) \, (9.69) 

for all k > 0. Hence by the same argument of non-stationary phase, we obtain the 
bound t _1 (|t - r\ + 2~ 2 H)~ N for any N > 0. Finally, if \t - r\ ~ 2~ 2j t, then we 
integrate only once and in the region \p — p | > 5 G (0, 1), where ~ \t — r\ and 
~ |t-r|(p-p ), so 

r)| < 5 + jf J^jL dp < 5 + (|t - r\6y\ (9.70) 

Choosing 5 = \t — r]" 1 ^ 2 , we obtain the claimed estimate in this region. □ 

We can now state a key dispersive estimate. 
Lemma 9.6. Let 1 > v > and a > = + -. // *C 1, i/ien 



\\e a ^U(t,0)<p\\ L z L ,,-« < „ 
e^ f /(t, S )/( S )d S || L?LrCT < P7|| i? ^ CT+i ^ r / 2 - 



(9.71) 



Proof. The short-time interaction part is estimated simply by the conservation 
and integration in t: 

|/f (/,<?) < M\\f\\ Ll J\g\\ Llx , (9.72) 
Hence for any 1 > v > 0, we choose M = 2^ and a as above, whence 

2 2 ^M^ a < M = 2 ju (9.73) 
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Thus Lemma 19.51 implies that 

\I^9)\<r^\f\\ Ll J\g\\ Llx . (9.74) 

Using (I9.55p . we obtain 

| jj o (e^-^U(t jS )wf(s)\w9(t))dsdt\<^2\I^f,g)\ 



(9.75) 



3=0 

oo 

< £ £ 2-||A fc /|| L? J|A„|| L ^ + ^-nfhijgh 

3=0 |i-fc|+|i-^|<2 

< ll/llifH-llfl'llLfia, 

which implies through duality 

' ^•»*U(t,s)f(8)d8 < \\V v {xYf\\^ m ~ PVIIxfLl^ ( 9 - 76 ) 

L t L x 

and the same estimate for the integral on t < s < oo. The remaining estimates 
follow by the usual duality argument. □ 

Lemma 9.7. For o > 14, one has 

!l(r -^|| ifi w<||a||l r |b|| i?4 ,- CT . (9.77) 

Moreover, if \\a\\L™ is small enough, then 1 — T is invertible on LfLl , ~ cr as well as 
on LfV u ^L 2 x ~ a forv>0 small. 

Proof. Decompose T — T into the long and short time interactions as before: 
(T - T)g = T s g + T L g, 

Tsg :=E^ = E f e (*-^°(l - U(t, s^AjVgis) ds, (9 ' 78) 

where Mj ^> 1 will be chosen later and Aj is as above. T 3 L is defined in a similar 
fashion. We estimate T 3 S in two ways: on the one hand, by definition of U(t, s), and 
the fact that V gains a derivative, 

Wet-^il-U^s^VtpWv < \\{l-U{t,s))Vip\\ L * 
< \b(t.,s)\\\Wip\\ L 2 < \t - s|||a|| L ^||^|| L 2,- CT , 

for any a, and on the other hand, 

|| e (t-.)i£o(l _ u(t, sfiAjVipWv < 2|||Vr 1 ^ J -W^|| 2 < 2-i\\<p\\v,-. (9.80) 

Hence, 

|| e (i- s) ^ 0(1 _^ (M))A ^|| L2 < It-slh-^WaWl^WipW^,-, (9.81) 
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By integration in s, for any uGl, 

\\ng\\ LfL 2 < M/2-i||a||l r |b|| i?i ,,- CT . (9.82) 

For T° L we use Lemma 19.51 with and without a (or a = 0), and a large. Taking 
account of the derivative gain from V and its decay, we obtain 

ll^|| L? ^<2 2 ^ f ^lbll^-- ( 9 - 83 ) 

Hence, by choosing Mj = 2s\\a\\ L ^, a > 14, and summing in j > 0, we obtain 

\\( T s + T L )g\\ L 2 L 2,-a < \\a\\l ? \\g\\ qLl ,-«. (9.84) 

Since 1 — T is invertible on L\L 2 ^~ a by Lemma 19.41 so is 1 — T by the Neumann 
series, if IHIl 00 is small enough. □ 

Now fix v > small and a large. Going back to the equation (19.501) of z, we 
obtain via the previous lemmas, 

WD-^zW^,-. < \\V-»/ 2 z \\ L2L 2,-« + \\V^ 2 F\\ L2L ^ 1h - v/2 

(9 85) 

| L 2- D -„/2 L 2, CT+L i L 2, 



< IU(0)|| 2 + \\a\\ L ~\\V~" /2 z\\ T 2 T 2,-* + ||F| 



and so, for ||a||L^ small enough, 

11^11^2^-/2^2,-^ < ||z(0)|| 2 + \\F\\ L iD-»/ai£.° + Lli£- (9.86) 



On the other hand, applying the Strichartz estimate for e lS,ot (see, for example, [191 
Section 4]) to (19~48|) . yields 

\\A\stz < \WWstz + \\Vz + F\\ Stz * < \\z(0)\\ 2 + \\z\\ L2v „ /2L 2,-<r + \\F\\ S tz*, (9.87) 

t x 

where the Strichartz space is defined by 

Stz := Lfh\ n L 2 B^ /6 . (9.88) 
Thus, we obtain the following estimate as a conclusion to this section: 

\\z\\stz < \\z(0)h + \\H L i Ll+ q {B l/l^ v -^ L 2,. )x - (9.89) 
As far as scattering is concerned, we note that 

poo 

||« — ^oolla — >• (t^oo), z oo = z + e i{t - s)2o U(t,s)F(s)ds. (9.90) 

Jo 

Hence there is ip + G L 2 such that 

||f/(0,t)z(t)-e i£o V+l| 2 ^0. (9.91) 
The following result summarizes the dispersive estimates of this section. 
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Proposition 9.8. Let v > be small and a large. Then there exists a small 5 > 
such that if < 5, then the solution of 

lt = iV/:j + P c [a{t)-Vj + f], T (0)=P c7 (0), (9.92) 
satisfies on any interval I 3 0, 

II^IL^LlnL^B^/W^L 2 .--)* ~ ll"K°)l|2 + \\f\\LlLi+L 2 t {B^l 2 r\V-» / 2 L^) X > (9.93) 

and for some 7+ G L 2 , and b(t, 0) = f* a(s) ds, 

\\j(t,x-b(t,0)) -e im 7+l|2 ->0 (f->oo). (9.94) 
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